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We develop a nonlinear generalisation of the causal linear thermodynamics of bulk viscosity, 
incorporating positivity of the entropy production rate and the effective specific entropy. The 
' \£} ' theory is applied to viscous fluid inflation (which is necessarily far from equilibrium), and we find 

, thermodynamically consistent inflationary solutions, both exponential and power-law. 
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£h ! T. INTRODUCTION 

in _ 

Scalar dissipative processes in cosmology may be treated via the relativistic theory of bulk viscosity (see |lj , 0] and 
references cited there). The causal and stable thermodynamics of Israel/ Stewart provides a satisfactory replacement 
of the unstable and non-causal theories of Eckart and Landau/ Lifshitz. However, it shares with these theories the 
limitation that it is based on assuming small departures from equilibrium, so that the transport equation is linear 
in the bulk viscous pressure. While this linear assumption is reasonable for many cosmological and astrophysical 
£N) ■ situations, there may be processes which involve large deviations from equilibrium. Such deviations are likely to lead 
to a breakdown of the theory. (This is known to occur for heat flux pL Q].) 

For example, if viscosity-driven inflation occurs (leaving aside for the moment various questions about the hy- 
drodynamical consistency of such models), then this necessarily involves nonlinear bulk viscous pressure Q]. Any 
application of the linear theories to viscous inflation (see Q for a review), unavoidably requires an assumption, which 
• is usually not made explicit, that linear thermodynamics remains valid well beyond the regime in which it is derived. 
Oh| The results of such a linear analysis are questionable. The alternative is to develop a nonlinear generalisation of the 
Israel/ Stewart theory for a more satisfactory model of viscosity-driven inflation. The goal of studying viscosity- 
driven inflation is to find inflationary solutions arising strictly from 'ordinary' matter in a far-from-equilibrium state, 
without invoking scalar fields. Of course, such models also need to address the question of how the observed density 
perturbation spectrum can be produced if one eliminates the usual scalr inflaton. This important problem is not 
\ addressed here. 

We present a phenomenological model of nonlinear bulk viscosity and apply it to viscous inflation. In the absence 
of a suitable microscopic foundation for a nonlinear theory, we propose a model which at least satisfies reasonable 
physical constraints, in that it: 

_ i 

1. approaches the Israel/ Stewart theory in the linear regime (i.e. small deviations from equilibrium), thus ensuring 
causal and stable behaviour in that regime; 

2. ensures that the entropy production is non-negative, i.e. the second law of thermodynamics is satisfied; 

3. imposes an upper limit on the bulk viscous stress, analogous to existing nonlinear generalisations of the heat 
flow equation |5j . 

In Section II we discuss briefly the problem of expanding fluids that are far from equilibrium, especially the case 
of inflationary expansion. On the basis of the above guiding principles, we then put forward in Section II the 
mathematically simplest form for a nonlinear transport equation, and consider some of its properties. Hopefully this 
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simple model can illuminate some of the overall features of a microscopically-based nonlinear theory, which would 
probably depend heavily on the particular interactions involved. We find in Section III exact solutions in a flat FRW 
universe, for simple and consistent thermodynamic coefficients and equations of state. Provided the thermodynamic 
parameters satisfy certain conditions, the model admits thermodynamically consistent inflationary solutions. Although 
viscous fluid inflation is the only application considered here, the nonlinear theory could in principle be applied to 
other far-from-equilibrium processes. 



II. TRANSPORT EQUATION FOR BULK VISCOSITY 

The hydrodynamic description which is implicit in fluid dynamics, requires that the mean interaction time t c of 
fluid particles should be much less than any characteristic macroscopic time-scale. In cosmology, this requires 

t c « H- 1 , (1) 

where H is the Hubble expansion rate. Condition ([!]) means that the fluid is able to adjust to the cooling caused by 
expansion, and to establish a state amenable to a hydro-thermodynamical description. In particular, if the fluid is in 
or close to equilibrium, then (Q) should ensure that the fluid has a well-defined local temperature. 

For a fluid far from equilibrium, there will not in general be a well-defined temperature. For example, very high 
and rapidly varying temperature gradients could give rise to large heat flux, without a meaningful temperature at 
any event. This could happen without violating (yj). 

If we focus on scalar dissipative effects, in particular ruling out heat flux, then large deviations from equilibrium 
arise from large bulk viscous stresses, i.e. |IT| > p, where p is the local equilibrium pressure, and the effective 
non-equilibrium pressure is 

p eff =p + IL. (2) 
In an expanding fluid, the dissipation due to II leads to a decrease in kinetic energy and therefore in pressure, so that 

n<o. 

Qualitatively, bulk viscosity is the macroscopic expression of microscopic 'frictional' effects that arise in mixtures 
H , |7j . The bulk viscosity tends to be largest when the contrast between components of the mixture (specifically, the 
contrast in cooling rates) is greatest. For example: 

1. a mixture of ultra-relativistic and non-relativistic particles in an ideal gas in the intermediate regime mc 2 ~ kT; 

2. a mixture of photons with long mean free path, and thermalised non-relativistic electrons and protons, whose 
mean free path is very short (radiative hydrodynamics); 

3. a mixture of gauge bosons which have acquired mass after a phase transition in the very early universe, and 
massless and effectively massless thermalised species. 

Bulk viscosity allows us to describe a mixture of different species effectively as a single fluid, provided a hydrody- 
namic description is reasonable. For example, in radiative hydrodynamics, if the thermalised matter dominates the 
energy density, then since it obeys ([!]), we expect that a single-fluid model is reasonable, even though the photon 
component strongly violates (Q). Far from equilibrium, the fluid model itself may break down, in the sense that (|l|) 
may not be satisfied by any interactions. We will assume that the deviations from equilibrium are not such as to 
cause a breakdown in the fluid description. The non-equilibrium conditions could lead to a weakening of ([!]), i.e. to 
t c <H-\ 

In the case of viscous fluid inflation (discussed in more detail in the following section), it is difficult to see how 
short-range interactions in the fluid could maintain a rate that is greater than the inflationary expansion rate. Fluid 
particles are being separated at tremendous speed, so that ([!]) should be strongly violated and a hydrodynamic 
description should break down. Any consistent fluid model of inflation needs to provide a way around this problem. 
One possibility could be the existence of long-range interactions, possibly involving a coupling of particle interactions 
to gravity, whose rate remains higher than the inflation rate. In our subsequent discussion of viscous fluid inflation, 
we will assume that some such mechanism exists to ensure consistency of the model. (Without such a mechanism 
being explicit, the model remains phenomenological.) Our task is then to provide a nonlinear transport equation to 
deal with the large bulk viscous stress involved in inflation. 

A covariant approach to the causal thermodynamics of relativistic fluids |9j is based on the hydrodynamic 
tensors n a (particle number 4-current), T a p (energy-momentum tensor) and S a (entropy 4-current), which are 
subject respectively to number and energy-momentum conservation 
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n% = 0, T Q/3 ;/3 = 0, (3) 

and to non-negative entropy production rate (the second law of thermodynamics) 

S%>0. (4) 

The entropy 4-current is assumed to be algebraically determined by n a and T a p alone: 

S a = S a (n f3 ,T^). (5) 

This is a reasonable assumption near equilibrium, and is well-motivated by kinetic theory B]. Far from equilibrium, 
(|J) could be a drastic assumption, in the sense that the full structure of non-equilibrium states is unlikely to be 
captured by only the first two moments of the particle distribution. However, for only scalar dissipation, it may be a 
less restrictive assumption. Furthermore, in homogeneous FRW spacetimes, the assumption that S a does not depend 
on spatial gradients of the moments will not be restrictive. In any case, without a clear alternative, we can begin by 
assuming that (||) is also valid for nonlinear scalar deviations from equilibrium. 

At each event, one can define a local reference equilibrium state, characterised by number density n, energy density 
p, pressure p, specific entropy S and 4-velocity u a . The arbitrariness in the reference state can be used to match n 
and p to the actual state, in the particle frame, which is characterised by vanishing particle flux, so that 

n a = nu a . (6) 

The energy-momentum tensor for scalar dissipation is then 

T af3 = pu a up + {p + II)/i Q| a , (7) 

where h a p = g a /3 + u a u/3 is the projection tensor, with g a p the metric. In the absence of vector and tensor dissipation, 
the entropy 4-current will be of the form 

S a = S eff n a , (8) 

where S e ff is the effective, non-equilibrium specific entropy. 
The conservation equations (||) for (||) and (Q) are 

n + 3Hn = , (9) 
p + 3H(p + p + n) = 0, (10) 

(p+p + n)u a + V(p + n),/9 = o, (ii) 

where 3H = u a ia is the average volume rate of expansion (H = Hubble rate in FRW spacetime), and ii a = u a - t pvP is 
the 4-acceleration. 

The transport equation for n arises from imposing the second law (|j) on (||) - which requires us to specify the form 
of S e ff. In the Israel/ Stewart theory 

%/ = *-(^)n 2 , (12) 

where n) is the bulk viscosity and r(p, n) is the characteristic time for linear relaxational or transient effects. The 
local equilibrium variables S and T satisfy the Gibbs equation 

TdS = (p + p)d( - ) +~dp, (13) 
V n I n 



which implies, using (|9j) and (10), that 



It is a feature of the Israel/ Stewart theory that dissipative contributions to the entropy 4-current are quadratic. Any 
linear contribution is ruled out by the requirement that the entropy density be locally a maximum in equilibrium. 
In the Israel/ Stewart theory, the deviations from equilibrium are small, i.e. 



3 



|n|« P . (15) 

Since nT ~ p and t/( ~ p -1 ~ p -1 , it follows from ( p"2| ) and ( |i~5| ) that 5 — S^// < 1, so that (|l5|) prevents the 
effective specific entropy from becoming negative. We will not impose the constraint ( |l5| ) on |II|. However, for 
simplicity, we will retain the quadratic form ([l2]), i.e. we will not add cubic and higher contributions, since there is no 
clear indication as to what such contributions would be - and since they would introduce further, unknown coupling 
coefficients. Since we require non-negative S e ff, a consequence of adopting ( |l2] ) in the nonlinear case is that the bulk 
stress is bounded: 



s ef f>o «. |n|<|n|_ = ^^. (16) 

Apart from this new feature, the main difference between the nonlinear theory presented here, and the linear Israel/ 
Stewart theory, reduces to the different ways in which the second law is imposed. 
By (|), © and (gjj), we find that (|) and (£|) lead to 

S a ;a = --X where 

r • t ( t C f\ 

X = 3H+-U+-U[ 3H+--±--\ . (17) 

In the Israel/ Stewart theory, the second law (Q) is satisfied by assuming a linear relation between the thermodynamic 
'force' X and the thermodynamic 'flux' II, i.e. 

n = -(X (Israel/ Stewart) . (18) 



(The linearity may be justified via near-equilibrium kinetic theory |8J.) Formally, it follows from ( 18[ ) that |II| is 
unbounded in response to growth in X. The near-equilibrium condition (|l~5| ) has to be taken as an additional, 
extraneous constraint on the transport equation (|i~8|). This is a sign of the limitations of the theory arising from 
its linearity. Beyond near-equilibrium conditions, the linear transport equation (|l^) cannot be expected to give 
reasonable predictions. 

The nonlinear generalisation of (|i~8| ) which we propose must be consistent with the upper bound (|l6| ) on |n|, so 
that while |n|/p is not restricted to be small, it does reach a limit as the thermodynamic force is increased without 
bound. This seems to be physically reasonable, and is in line with existing phenomenological generalisations of the 
heat transport equation The simplest generalisation of ( |i~8| ) that we could find with this property, and the other 
properties listed in the introduction, is 

CX 

II = (nonlinear) , (19) 

1 + T*X 

where t» > is a characteristic time-scale for nonlinear effects. Nonlinear effects are significant for X > t~ , 

For small thermodynamic force X, equivalently small |n|, ( |l~9| ) tends to the form (|l^), and the linear Israel/ Stewart 
theory is recovered. For \X\ — ► oo, —IT tends to the upper bound C/ T *- It follows from ( p~7| ) and ( |i~9| ) that 



S a . a = nS e ff - — 



n 2 r 



1 + ? n 



(20) 



This means that the second law holds identically by virtue of the upper bound on the bulk stress: 



s a . a >o ^ -n<— . (21) 



As — II approaches C/ T *> the entropy production rate ( p0| ) tends to infinity, providing a barrier against values —II > 
(yV*- Note also that ( p0| ) implies that the rate of growth of effective specific entropy is greater in the nonlinear than 
in the linear theory. The condition ( p6| ) implies by (|l]) that 



T > . (22) 

* " v 2nTS V ' 



By (JT^) , the detailed form of the transport equation ( |l^ ) is 
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= -3(H- ±tU 



II(1 + 3t*#) 



t C T 



1 + ? n 



(23) 



The Israel/ Stewart theory, in its full (non-truncated) form is the limiting case r* = 0, when the terms in 

round brackets in ( |2^ ) all reduce to 1. 

The thermodynamic coefficients C, and r are known from the linear theory, based on kinetic theory arguments Q . 
For the nonlinear coefficient r» we do not have a corresponding microscopic derivation. However the phenomcnological 
theory presented here imposes on r* the limit (22), which is determined by the linear coefficients and local equilibrium 
variables. 



III. THERMODYNAMIC ALLY CONSISTENT INFLATION 



We specialise now to a flat FRW universe 

ds 2 = -dt 2 + a(t) 2 (dx 2 + dy 2 + dz 2 ) 



where H = a/ a. The transport equation (|23|) implies, with the Einstein field equations, an evolution equation for H, 
once we specify the thermodynamic coefficients and equations of state. The field equations are 

H 2 = \p, (24) 
H + H 2 = -±(p + 3p + 3II) . (25) 

For the equilibrium pressure we impose the linear barotropic equation of state 

p=(y-l)p, (26) 

where 7 is a constant, I < 7 < 2. If T is also barotropic, T = T(p), then the integrability condition d 2 S/dpdn = 
d 2 S/dndp of the Gibbs equation (||) leads to @ 

T oc p^-^h (27) 

We note that the ideal gas law p = nT is incompatible with barotropic temperature and the linear barotropic equation 
of state (^), unless II = jo]. This follows from ( p7| ) - which is a consequence of ( p6| ) and T = T{p) - and the 
conservation equations (^) and (|To)) . 

The linear relaxation time is related to the bulk viscosity by |To| 

" " (28) 



v 2 (p + p) v 2 "fp 

where v is the dissipative contribution to the speed of sound V , so that V 2 — c 2 +v 2 , with c s the adiabatic contribution. 
By causality, V < 1, so that 

v 2 <l-c 2 s =2~ 1 , (29) 

where we used (p6|). We assume that v is constant, like c s is when ( |26| ) holds. The bulk viscosity itself is often taken 
to be of the simplified form ( oc p q l 2 (q constant), which should not be unreasonable in FRW spacetime with j2^). 
By (24) we have 

C = aH q , (30) 

where a (> 0) is constant. Finally, we need to specify the nonlinear characteristic time r + . The simplest possibility 
appears to be 

T*=k 2 T, (31) 



where k is constant, and k — is the linear (Israel/ Stewart) case. This ansatz is subject to the condition (|22j), which 
does not take a simple form in general. 
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From ( |24| ) - (|26|) we get 

n = -2H - 3-/H 2 
Together with @ - (|l]), this brings (§|) to the form 



(32) 



> >2 



/ 2fc 2 
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H + 3HH 



3jv 2 
~2a 



ak 

"fV 2 



Hi- 



ll 
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(33) 



This is the fundamental dynamical equation for nonlinear bulk viscosity in a flat universe. Once (|33|) is solved for H, 
p is determined by (Q) and II by |3^). When k — 0, the terms in square brackets all reduce to 1, and we recover the 
dynamical equation for the (non-truncated) Israel/ Stewart theory jl| .[] 

Viscous inflation has been discussed in a number of papers (sec |l| for a review), both as a phenomenological model 
of quantum particle and string creation effects, and as a fluid alternative to scalar-field inflation. Particle production 
is formally equivalent to an effective bulk viscosity ||, but we will confine ourselves here to discussing fluid models, 
in which the bulk viscous stress | H] has become large enough to make the effective pressure (Q) negative and initiate 
inflation, and where net particle creation and decay occur only during reheating, at the end of inflation. (Causal linear 
thermodynamics, incorporating net particle production, may be used as a phenomenological model of reheating p2[.) 
As discussed in the previous section, we assume that long-range interactions or some other microscopic mechanism 
are present to ensure that a hydrodynamical description remains reasonable during inflation. 

The condition for inflation a > implies by (p5|) that 



n > P 



(34) 



It is clear from ( [34] ) that viscous fluid inflation is not close to equilibrium, i.e. it strongly violates the condition 
(H), and we need to use the generalised nonlinear equation ([53]). We consider first whether ( |33| ) admits exponential 
inflation, and if so, whether the solutions conform to the thermodynamic conditions discussed above. 
Setting H = H , (|35| ) gives 



H, 
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277/ 



(2v 2 - 1) (v 2 - k 2 ) 



By (32) and (E4h, n = -37^ = -7/). Then the condition (21) implies, using <M) and (RQ), a limit on k: 



k < v . 

In order to check the condition (|2|), we evaluate the specific entropy. Integrating (jlj) with T = Tq = constant: 



S = S + 



V n T 



where n = n^a . Then (12) leads to 



Seff = So 



37^ 

2n T v 2 



(2v 2 - l) e 3Hot . 



Thus we require v 2 > h. It follows from (p5|), ( pq ) and ( P9| ) that exponential inflation is possible provided 

k < v , i<i; 2 <2-7, 



(35) 



(36) 



(37) 



(38) 



(39) 



and then the entropy generation and effective specific entropy are positive. Thus the solutions are thermodynamically 
consistent. Note that 7 < I. Exponential growth of S is also found in the linear theory fl - but the advantage here is 



1 Note that [jjj effectively assumes v 2 = I/7, and takes the exponent in (^) to be a free parameter r; one must set r = (7 — 1)/7 
10- 
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that the theory is nonlinear, and we consider the effective specific entropy, which is more suitable far from equilibrium. 
The apparent closeness of the linear and nonlinear results is an artefact of the constancy of H and thereby of p, T, II. 
As shown in pj], it is possible to generate a large amount of entropy via the exponential form (|38|). 

Following the qualitative analysis techniques for the Israel/ Stewart solutions jn]], jyj), we find that the stability 
properties of the solution (pq) depend on the parameter 



2v 4 



2k 2 



4 7 ( v 2 - k 2 ) z (2v 2 - 1) 



> 0. 



(40) 



Then the fixed point (0, Hq) in the phase plane {(H, H)} has the following properties: 
for 1 — u> < q < 1, it is an asymptotically stable node; 
for q > 1, it is a saddle point; 
for q < 1 — oj, it is an asymptotically stable focus. 
Thus the solution is an attractor for q < 1. 

When q — 1, exponential inflation can occur with Hq arbitrary, provided that the bulk viscous parameter satisfies 



2 7 v 4 



{2v 2 - 1) {v 2 - k 2 ) ' 

as follows from ([55]). Of more interest is the existence when q = 1 of power-law solutions a 
inflationary for N > 1. By (|33"1), these solutions exist if 

277 3 [a (1 - 2t> 2 ) (v 2 - k 2 ) + 2 7 i; 4 ] TV 3 
+ I87 2 [3ak 2 - 2 7 i; 4 - 2av 2 (l + fc 2 )] iV 2 
+ 12a 7 (v 2 - 3k 2 ) N + 8ak 2 = . 



a t 



(41) 

1 which are 



(42) 



Clearly there are solutions of this cubic for a wide range of thermodynamic parameters. However, we are concerned 
to identify the thermodynamically consistent solutions, for which the conditions are stringent. 

By (jH), n = -N(3jN -2)t~ 2 , and together with (§J), @, (§§) and @, this shows that the condition @ for 
positive entropy generation is 



Integrating dll) we get 



and then (O) gives 



fc < 



5 — 50 



3 7 7V 
37iV- 



1/2 



(43) 



,1/7 



7(3A^ S 



^(37^-2)77 



5o 



37V^ 



,(l-7)/7 



{9 7 2 (2« 2 - 1) N 2 + 12jN - 4} t( 3 ^-2)/7 , 



The effective specific entropy is positive if the ./V-quadratic in braces is positive. This requires that 

1 



v < 



V2' 



and it places an upper limit on N, so that 



2 2 

— < N < ; 

37 3 7 (1 - 



V2v) 



(44) 



(45) 



(46) 



(The lower limit ensures that II < 0.) 

Thus the constraints ((43|), ( (45| ) and ( p6| ) must be satisfied for any solution of ([l2]) to be thermodynamically consistent. 
We present a simple example to show that consistent inflationary solutions do exist: 



7= 3' 



= \=k\ a = %, N = 2. 
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In conclusion, we have argued for a nonlinear generalisation of causal linear thermodynamics in order to describe 
viscous inflation (without particle production). A simple phenomenological theory, supplemented by a simple model 
for the thermodynamic coefficients and equations of state, leads to severe restrictions if thermodynamical consistency is 
imposed. However, consistent solutions do exist, both exponential and power-law. These solutions respect the second 
law of thermodynamics and have positive effective specific entropy. Of course, the phenomenological theory lacks a 
microscopic foundation, and it is necessary to assume a priori that long-range interactions or some other mechanism 
are present in order to maintain the viability of a thermo-hydrodynamical description during inflation. If the theory 
presented here were applied to non-inflationary processes that are far from equilibrium (e.g. in astrophysics), then 
this limitation would not apply. 
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